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Thermodynamics of a heavy quark-antiquark pair in SU(3)-QCD is studied both 

■ below and above the deconfinement critical temperature T c . In the quenched case, 

O , 

a model of the string passing through heavy valence gluons yields a correct estimate 
H of T c and the critical behavior of the string tension below T c . For two light flavors, 

entropy and internal energy below T c can be obtained from the partition function of 
heavy-light mesons and baryons. To calculate the free energy of the system above 
T c , we apply second-order perturbation theory in the interaction of the quark-gluon 
\^ \ plasma constituents with the static quark-antiquark pair. The results for the entropy 

in ' 

CNJ ■ and internal energy, obtained both below and above T c , are compared with recent 
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lattice data. 

I. INTRODUCTION 

Recent RHIC experiments suggest that the quark-gluon plasma may be more 

like a perfect liquid [5|, where the mean free path of a particle is much smaller than the inter- 
particle distance. Values of the mean free path following from perturbative calculations jf|, 
however, are large. This makes perturbative approaches to the plasma questionable. 

Lattice simulations of the quark-gluon plasma also need non-perturbative methods for 
their explanation. For instance, it has recently been argued [7fl that non-perturbative 
chromo- electric fields can describe the singlet free energy of a static quark-antiquark pair at 
large distances which remains quite sizeable above the deconfinement temperature T c . Other 
lattice data calling for a theoretical explanation are the anomalously large maxima of the 
internal energy and entropy of the static quark-antiquark pair in unquenched QCD around 
T c 8] (see Ref. [s| for a review). Equally surprising is the rapid fall-off of these quantities 
right after the phase transition. This paper aims to analyze what kind of nonperturbative 
physics is necessary to understand these lattice simulations. We consider three cases: In 



finite-temperature quenched QCD the QQ system contains many valence gluons 
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which may form a gluon chain. We demonstrate that the high entropy generated by this 
object gives a reasonable prediction for the critical behavior of the string tension near T c . 
In unquenched QCD for T < T c the QCD string may break and the produced light quark 
and antiquark interact with the heavy QQ pair to form mesons. In unquenched QCD for 
T > T c the constituents of the quark-gluon plasma interact via screened gluon exchange with 
the static quark and antiquark. These three simple models are supposed to catch the main 
physical aspects of the heavy-quark system in hot QCD. They cannot represent a unified 
picture due to their inherent approximations. By comparison with the lattice simulations 
one sees how realistic the so obtained picture is. 

The paper is organized as follows. In the next section, we consider quenched SU(iV c ) 
QCD below T c , calculate the partition function of a gluon chain, and make predictions for 
T c and the effective string tension cr(T). In section III, using cx(T), we calculate within 
the relativistic quark model the partition function of heavy-light mesons and baryons. The 
entropy and internal energy stemming from this partition function are compared with the 
corresponding lattice data. In section IV, we calculate the interaction energy of quarks, 
antiquarks and gluons with the heavy quark-antiquark pair immersed in the plasma. We 
impose that the plasma particle-number densities vanish at T = T c . With this constraint 
implemented we calculate the change of the entropy and internal energy in second order 
perturbation theory due to the quark and antiquark pair. Finally, we estimate contributions 
coming from possible light qq octet bound states, which may exist after deconfinement. In 
section V, we summarize the main results of the paper. 



II. STATIC QQ-PAIR WITH GLUONS AT T <T C 

The QQ-string which sweeps out the flat surface of the corresponding Wilson loop is 
produced by soft stochastic gluonic fields. In addition, fluctuations of the gauge field exist, 
which lead to string vibrations. These fluctuations may be related to valence gluons through 
which the QQ-string passes. At asymptotically large QQ-separations considered here the 



string passes through many valence gluons and forms a gluon chain 12]. Further, the 
energy of one string bit between two nearest gluons in the chain is constant. As long as 
the thermal mass of a valence gluon is smaller than this energy, gluons move together with 
the string and do not affect the global dynamics of the string (see Fig. 1). However, at a 
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Figure 1: Gluon chain for T < Tq where valence gluons move together with the string. 








Figure 2: Gluon chain for T > Tq where the valence gluons become static, but color may change 
from one string bit to another. 

certain temperature T the gluon's thermal mass (oc T) becomes larger than the energy of 
one string bit. From this temperature on, the system looks totally different, since gluons 
from the string's standpoint are now nearly static. Therefore at T < T < T c a gluon 
chain becomes a sequence of static nodes with adjoint charges, connected by independently 
fluctuating string bits (see Fig. 2). At the moment of formation of such a chain, its end-point 
originating from the heavy Q performs a random walk towards Q over the lattice of static 
nodes. The entropy of such a random walk turns out to be large, namely proportional to 
its length, and eventually leads to the deconfmement phase transition in this model. The 
reason for a large entropy is that color may alter from one node of the random walk to 
another, i.e. every string bit may transport each of the N c colors. A similar picture has 



been proposed by Greensite and Thorn 



12| as a time cut of a large Wilson loop leading to 



the above sequence of valence gluons arising from double lines of neighboring plaquettes. 
Therefore, the total number of states of the gluon chain grows exponentially with its length, 
L, as N^ a , where a is the length of one bit of the chain. 
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To quantify these considerations, let us begin with considering a free random walk, which 
is performed over a <i-dimensional Euclidean hypercubic lattice (in our case d — 4) with the 
spacing a. At zero temperature, a walker, starting from the origin, arrives at a distance R 
(R = \1Z\, where TZ is a d-dimensional vector) with the probability [13] 

where s is the proper time of the walk, i.e. s = a 2 x (number of steps). The proba- 
bility obviously obeys the conservation law / d d lZP(s, R) = 1 and the initial condition 

lim P(s,R) = 8^ (TV). This probability is related to the proper-time representation of the 

r(-— i] 

Green function (-<9 2 )^) = 4 J/l Rd L 2 through 



Let us generalize these formulae to temperatures below the temperature of dimensional 
reduction. The probability then depends also on the number n of a Matsubara mode. In the 
same way as at zero temperature, it can be obtained from the proper-time representation of 
the Green function (— d 2 )^o- Instead of Eq. j2j), one has (fixing from now on d — 4): 



CO 



-d 2 )n,o = £ / dsP n = £ /°° ^ exp 



R 2 + (Kt - (5n) 
4s 



(3) 



/O (47TS) 2 

where n in runs from — oo to +oo, and j3 = 1/T. Note that this representation can be 

n 

obtained directly from upon the decomposition R 2 — >• R 2 + (7£ 4 — (3n) 2 . In the form Q, 
the probability obviously obeys the conditions 

/ d 4 TZPJs,Tl) = 1, \imP n (s,TV) = 5 {3) (~R)5(K 4 - pn), 

which are similar to those at zero temperature. 

In case of a heavy quark-antiquark system, the random walk is not free, as "the walker" is 
attached to the heavy quark Q by the confining string. The length of the string is L = s/a, 
i.e. the number of steps s/a 2 times the step size a. Accordingly, Eq. ((31) becomes replaced 
by the Green function of the confined walker: 

ds 



G(R,T) = £/ dsP n (s,n)e-^ a = J2 7T-^exp 
„ Jo „ Jo (Airs) z 



R 2 + (pn) 2 pas' 
r — • 4 

4s a 



Here we use Eq. ([3D with 7£ 4 = 0, since the quark and the antiquark at the ends of the 
string are infinitely heavy. We will also use the zero-temperature value of the string tension 
a = (440 Me V) 2 . 
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The full partition function of the random walk differs from this expression by the above- 
mentioned entropy factor, N^ /a , which should be included in the s-integration. The expres- 
sion for the partition function is thus 



2(R,T) = W C 
™ Jo 



ds 



(Airs) 2 

Performing the s-integration, we obtain 



exp 



R 2 + (Pn) 
4s 



ap 



In N r , 



(5) 



Z(R,T) 



a 



\nN r 



4tt 2 \ a \T 



E 



lri7V c 



(R 2 + (Pn) 



R 2 + (Pn) 2 



where K\ is a Macdonald function. Next, to get the static QQ-potential, one should as usual 
take the limit of asymptotically large QQ-separations, |R| — ► oo. There, the Macdonald 
function falls off exponentially, for which reason the sum can be approximated by its zeroth 
term alone. Furthermore, since the random walk is suppressed at T < T , its free energy 
vanishes at T = T , and therefore Z(R, T) should be normalized by the condition Z(R, T ) = 
1. The effective string tension is then defined through the full free energy of the system, 
which is the sum of the usual linear potential and the normalized free energy of the random 
walk: 



a(T) 



T Z(R,T) 
o — — In 



R Z(R,T ) 



R— >oo 




(6) 



An estimate for T c now follows from the condition that the argument of the first square root 
vanishes: 

Equating T c to the modern N c = 3 lattice value [8|, [9||, 270 MeV, we obtain for the effective 
length of one string bit a ~ 0.31 fm. This is larger than the minimal possible value of this 
quantity, a = 0.22 fm - the so-called vacuum correlation length jl^J], which defines the onset 
of a string-bit formation. As for the temperature T , it can be defined from the condition 
cx(T c ) = 0, which yields 

T = : T f ~ 130 MeV. 



In N c + 1 

An important finding of our model is the behavior 



a(T) 



T c -T at T^T C . 



(8) 
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Figure 3: The ratio a(T)/a at N c = 2 according to Eq. JBJ) (full drawn curve) and Ref. [17j (dashed 
curve) . 



The same critical behavior follows from the Nambu-Goto model for the two-point correlation 
function of Polyakov loops [l5|. 



Let us consider the limiting case when string bits cannot alter color, i.e. one should 
formally set in Eq. (El) N c — 1, that yields 



a(T) 



a 



laT 
a 



T 



The fundamental difference of this limiting case from the realistic one iV c = 3 is that here 
cr(T) ~ (T c — T) at T — > T c . The corresponding value of the critical exponent v = 1 defines 
the universality class of the two-dimensional Ising model, which by no means can be realized 
in four-dimensional quenched SU(iV c ) QCD. The same linear fall-off of <r(T) with (T c — T) 
one finds also in the Hagedorn phase transition and in the deconfinement scenario based 
on the condensation of lon g Cosed String s fl Q. It co m es as a m ete consequence of the 
formula cr(T) = o — ^ with the entropy S oc R. For long strings the entropy is just the 
logarithm of the number of possibilities to realize on a lattice with the spacing a a closed 
trajectory of length L. Specifically for a hypercubic lattice, the entropy is S = ~ -ln(2d— 1). 
Therefore, the free energy of such a closed string, F = oL — TS, vanishes at T c = ln ^d~i) ■ 
Numerically, in d = 4 dimensions, T c = 103 MeV, that is by a large factor of 2.6 smaller 
than the modern lattice value 270 MeV. The principal difference of our calculation is that 
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Figure 4: Entropy S(T) of the static QQ-pair in the SU(2) quenched QCD at R = 1.5 fm and 
temperatures close to T c ~ 304 MeV, calculated from Eq. j6|). 



we consider the random walk between two different points, taking into account the confining 
force between the final and the initial points of the walk along its trajectory. These are the 
two facts, which eventually lead to a different critical behavior. In conclusion of this section, 
we note that the idea to have a large entropy of the QQ-string in the vicinity of T c due to 
the valence gluons has recently been mentioned in Ref. {t|. However, the possibility that 
this mechanism leads to deconfinement has not been formulated in that paper. 

It should finally be mentioned that, in the string models discussed in the previous para- 
graph, as well as in the model developed in this paper, the free energy vanishes gradually 
at T — ► T c . Therefore, none of these models can describe first-order phase transition, which 
takes place in SU(iV c ) quenched QCD at N c > 3. They only can describe the phase transi- 



tion in the SU(2) quenched QCD, which was studied on the lattice in Ref. [17|. In that case, 
the mean-field critical exponent v = 1/2, found in this paper and in Ref. 1_5|, is changed 
to the 3d-Ising one, v = 0.63. This makes our values of cr(T) larger than those obtained 
in Ref. [17l |. Using, similarly to that paper, the value T c = ~ 304 MeV, we plot in 

Fig. [3] the ratio a{T)/a following from Eq. ((6]) &t^N c = 2. In the same figure, we present a 



curve interpolating the lattice data from Ref. [17J|. In Figs. [4] and EJ we plot entropy and 
internal energy of the static QQ-pair in the SU(2) quenched QCD at R = 1.5 fm, calculated 
from Eq. © by the formulae S(T) = -R^ 1 , U(T) = a(T)R + TS(T). These quantities 
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Figure 5: Internal energy U(T) in GeV of the static QQ-pair in the SU(2) quenched QCD at 
R = 1.5 fm and temperatures close to T c ~ 304MeV, calculated from Eq. (|6j). 



read 
S(T) 



aR 



P 



+ (3 C 



IP-Pc 



aa 



and U(T) = aR 



1 + 



2Jaa(P - P c 



2^/aa((3 - & 

respectively. In these equations, we use the value aa ~ 0.21 GeV, which stems from Eq. (|7j) 
upon the substitution N c = 2, T c = 304 MeV. This value of aa corresponds to a ~ 0.21 fm, 
which i s ag ain larger than the vacuum correlation length in the SU(2) quenched QCD, 
0.16 fm 18]. Both S(T) and U(T) have a singularity of the type (T c - T)- l l 2 at T -> T c . 
This singularity is weaker than 8 [ T 'j^ T ^j , which takes place for the first-order phase transition 
at N r > 3. 



III. STATIC QQ-PAIR WITH LIGHT QUARKS AT T <T C 

Let us now consider unquenched QCD with a heavy quark pair QQ at large separation 
(R > 1.5 fm). Then the QQ-string breaks due to the production of a light gg-pair. The 
subsequent hadronization process leads to the formation of a heavy-light meson (Qq) or 
a heavy-light-light baryon (Qqq), together with their antiparticles. We will calculate the 
thermodynamics of these hadronic objects as a function of temperature. We consider the 
(Nf = 2)-case, with light u- and <i-quarks, and use the value T c = 200 MeV given in in 
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order to compare with the corresponding unquenched lattice simulations. Equation ((71) then 
yields an effective value of a = 0.23 fm, which is now used in Eq. ([6]) to determine cr(T). For 
a heavy-light meson, the Hamiltonian of the relativistic qQ-system is 



h q, =™q+ vp 2 + m \ + v ^ with v ^ = a{ ?y - ( 2 - v a ( T )- (°) 

The parameter tuq denotes the heavy-antiquark mass, m q is the constituent mass of the light 



quark, m q ~ 300 MeV. The subtraction of2ya(T) in V(r) is known to be important to reach 
agreement between the predictions of the relativistic quark model with the phenomenology 
of meson spectroscopy 20, 21|. An additional correction 5^Ja(T) ensures that the mass of 
the lowest state in our calculation corresponds to the lightest heavy-quark meson. We use 
the D°-meson with m^o = 1.864 GeV and trq = m c = 1.48 GeV to fix the parameters of 
our model. 

The square root in the Boltzmann factor can be handled by integrating over an auxiliary 
parameter: 

exp (-/y P 2 + ml) = — J Q rf/xexp j^ 2 - — (p 2 + m\) j . 
One then arrives at a 3d Schrodinger equation of the form 

(-KV 2 + M|x|)V(x) = EiJ)(x), 
where K and M are positive constants of dimension [length] and [mass] 2 , respectively. Its 



eigenenergies read 



2ll 



E nr j = a nr>l (KM 2 ) 1 / 3 , 

where a nr /s are positive numbers with aoo — 2.34. In calculating the partition function of 
the heavy-light mesons (Qq and Qq) [2l| we take their ground state into account exactly 
and model their higher eigenenergies by those of a 3d harmonic oscillator with the frequency 

This strategy works well for the lowest states, and higher states are Boltzmann-suppressed, 
so the lack of accuracy does not matter so much. Because of the degeneracy factor (n/2 + 
l)(n + 1) for the oscillator, we have 



• / 2mes 



4Z? 
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with 

2 



Z Qq = -= exp {-(5m Q + (2 - 6)Pyfi(T?) x 



00 

X 

n=0 



C J_ ^ 2 P m q 

(n + a 00 )(Juj - fj, - -j-^- 



(10) 



The partition function of two noninteracting mesons is where the factor 4 is a product 

of two flavors and two spin states of the mesons because the spins and isospins of the 
produced q and q are coupled to isospin and spin by string breaking. The mass uiq of 
the heavy antiquark is dropped in Eq. (fTUl ) in order to compare our calculation with the 
lattice data jl^ . which simulate a heavy quark and antiquark by two Polyakov lines. Doing 
the sum over n analytically, we arrive at the following result: 
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F 2mes = —Tin (4 - 25) Ja(T) - 2T In 

7T v 



J dyuexp ( — ji 2 — 



P 2 m 2 \ exp(-a 00 f3u) 



A/j, 2 I (I- exp(-(3u)) 3 

(11) 

The remaining /^-integration is done numerically. The necessary value of the parameter 5 to 
reproduce the correct m D o mass is 5 = 0.354, which can be obtained from the free energy F 
in the limit T — > 0. Entropy and internal energy can be calculated by the standard formulae 
S = —dF/dT, U = F + TS. Note that, when calculating the entropy, the derivative d/dT 
should not act on the Hamiltonian in the partition function, since otherwise the consistency 
of the thermodynamic relations would be violated. 

It is also possible that string breaking generates a baryon and an antibaryon instead of 
two mesons. The Hamiltonian for a baryon reads 

2 

H Qqq = m Q + Y. (\/p" + m2 g + V{ri)) , 
i=i 

where we approximate the position of the baryon string junction by the position of the heavy 
quark Q, which is legitimate due to the heaviness of Q. Since the heavy quark Q is static, 
terms from the two light quarks separate and one can deal with this Hamiltonian similar to 
the meson case. The corresponding free energy of baryon and antibaryon reads 



F 2bar = -T In ^ - (8 - 45) J^T) - AT In 



TV 2 



J dfi exp ( —[i 



2 P 2m g \ exp(-a 0Q pu) 
Afi 2 J (1 -exp(-/3cj)) 3 

(12) 



For simplicity, we restrict ourselves to diquark quantum numbers with zero isospin and 
zero spin. The parameter 5 in the baryon calculation is adjusted so that the corresponding 
ground state corresponds to the lightest charmed baryon A+ with m(A+) = 2.286 GeV. The 
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Figure 6: The calculated entropy S(T) (full drawn curve) of two mesons and two baryons produced 
by string breaking for a heavy Q and Q at large separation (R > 1.5 fm) is shown as a function of 
T/T c with T c = 200 MeV. The stars show the lattice data [ijj] f° r the same QQ configuration. 




resulting value 5 = 0.393 is not very different from 5 = 0.354 obtained for the heavy-meson 
system, therefore we use the averaged value 5 = 0.37 in both cases. 

Combining the contributions to the partition function coming from mesons and baryons, 
we have 



In Figs. M and d we plot the entropy and the internal energy corresponding to Eq. ( fTBl 
together with the corresponding lattice data Here S(T) and U(T) mean the change 

in entropy and internal energy due to the presence of a heavy quark-antiquark pair in the 
hadronic heat bath. The heat bath justifies the calculation in the canonical ensemble in our 
case. The heavy Q and Q and the two light q and q from string breaking are the sources of 
our mesons. Therefore we think that a comparison of our model calculation with the lattice 
simulation is appropriate. The agreement is very good for the entropy, while the internal 
energy coming from our calculations is smaller than the lattice results. Our curve is shifted 
downwards by AU — —0.5 GeV compared with the lattice simulations. Due to the simplicity 
of our calculation the qualitative agreement is nevertheless good. 




F = -T\n(Z 2mes + Z 2bar ) = -Tin (AZ 2 Qq + Z± q ) . 



(13) 
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Figure 7: The calculated internal energy U(T) (full drawn curve) of two mesons and two baryons 
produced by string breaking for a heavy Q and Q at large separation (R > 1.5 fm) is shown as 
a function of T/T c with T c = 200 MeV. The stars show the lattice data [lj| for the same QQ 
configuration. 

IV. STATIC QQ-PAIR WITH LIGHT QUARKS AT T > T c 



In this section we consider a heavy QQ pair at temperatures T > T c , i.e. after deconfine- 
ment with free quarks and gluons present in the plasma. The heavy quark and antiquark 
are immersed into this plasma at the points xq and Xq at large separation (R > 1.5 fm). 
The non-Abelian interaction energy of a plasma constituent a with the QQ pair has the 
following form: 

V a (r) = C aQ U(\r - x g |) + C aQ U(\r - x |). (14) 

Here, C a ® and C aQ are the corresponding color factors, namely products of SU(3) generators, 
which are discussed in Appendix [Bj and the interaction potential U{r) is given by a screened 
gluon exchange between the colored sources 



W(r) 



4?T / ^ a *( k > T ) c *T m 



(2n) 3 k 2 + m 



(15) 



D 



For au , we use a running coupling constant at finite temperature T which has been derived 



in Ref. 



22, 3 



a s (k,T) 



k 



1 + exp (u 2 ^ - u 3 



13 



— ( ~ a " (16) 

l + exp(v 2 £ -v 3 ))ln(e+ (£) 1 / 



A, 



This running coupling increases for small momenta k, has a maximum at k ~ T and finally 
decreases for k — > oo. 
The Debye mass is 



This g is determined from the selfconsistency equation [221 ] with the running coupling. 

The interaction energy (TT4l) is to be summed over all the N constituents a of the plasma, 
V(r 1; . . . , r N ) = y^j^ 1 V ai (rA The perturbative expansion of the free energy in classical 
thermodynamics reads (see 24 ]) 

F = F + F QQ + (V(n, . . . , rjv )) - -L ((V 2 (r 1; . . . , rjv )) - (V(r 1; . . . , r^)) 2 ) , (18) 

where for each operator 



(0(n, . . . , rjv )) = Tr, ^ft / j 0(r!, • • • , r,v)e 



/3(F -E ) 



Here, Fqq is the free energy of the interacting heavy Q and Q, F Q is the free energy of the 
quark-gluon plasma, and the rest describes the interaction of the plasma constituents with 
the QQ-p&ir. Tri denotes the trace over color indices of the interacting plasma constituent 
in the QQ singlet state. The first-order term (V(ri, . . . , r N )) vanishes due to color neutrality 
of the plasma, the second-order term, however, yields a nonvanishing result and describes 
the various interactions of the plasma constituents with the QQ singlet state. The respective 
color factors can be calculated as products of generators in the various representations of 
SU(3). In second order, only the diagrams with two interactions of one plasma constituent 
with either Q or Q or both of them are nonzero. 



The energy of N q light quarks, Nq light antiquarks and N g gluons is E = \[ p 2 + mf, 



where N — N q + Nq + N g . Here the squared mass of quarks and antiquarks is 



25 



26] 



m 2 q = m 2 = ml + 2gT\\ ^1 m + gT^ ] (19) 



and the kinetic mass of gluons is given by m 2 = \m 2 D 



25 



at N c = 3. We also use g = 2.5 



as above and set the current quark mass to m = 30 MeV. 
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Although our heuristic introduction used classical statistics we will calculate in the follow- 
ing the free energy in quantum statistics as it is appropriate. The zeroth-order contribution 
F reads 



a=q,q,9 L ^ 

= -£iV a T In [AWTmlK, (^)) , (20) 



where K 2 is a Macdonald function. The contribution F Q of the quark-gluon plasma alone 
without the heavy QQ pair is subtracted in lattice calculations jj], Q. Therefore we will 
also leave it out in our results. 

The contribution of the QQ-interaction, 

4 e -m D R 

f qq( r , t ) = -gQ!« R » 

can be neglected because it is of the order of a few MeV at large distances R ~ 1.5 fm. 

The term (V 2 (ri, . . . , rjv)) in Eq. ( fT8l) is associated with two interactions of each plasma 
constituent. The contribution is the same for constituents of the same kind: 

F 2 = -^(V 2 ( ri ,..., rAr )) 

= -— (l n eS + -n cS + -n cS ) I d 3 a ^lH. (i _ ^(xq-xq)) (21) 
ttT U * + 9 q + 2 9 ) J % 2 + mlf V 1 6 ) W 

(see Appendix O). In this equation, the color factors c q = Cq = 2/9 and c g = 1/2 appear 

multiplying the effective densities rf^ = Nff /V of quarks, antiquarks or gluons, respectively. 

Since the heavy QQ pair is color neutral, the contribution F 2 vanishes, when the QQ - 

distance |xq — xq| approaches zero. The running of the coupling with temperature is crucial 

for the strong decrease of the entropy and internal energy with temperature. For a fixed 

value of a s we would have 

F ' = -Ml< + r? + l<) a ^( 1 - e ~ mDR ^ ^=|x Q -x | = 1.5fm, (22) 
and the corresponding entropy and internal energy would decrease with temperature even 
slower. 

The next step is to calculate the effective densities, starting from the free ones. The 
particle densities of free quarks and gluons are listed in Appendix 0, and the results are 

n q = n g = 3TN f — £ - K 2 i^—n 
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Figure 8: The calculated ratio P e s(T)/T 4 of Eg. ([281) (full drawn curve) is shown in comparison 
with the one from lattice simulations of Ref. [27y (dashed curve) for two light flavors. 



,2 oo 

12T— 9 - 



n=l 



-n 



(23) 



The full density is a sum of the three contributions n q + n q + n g . Confinement requires 



n 



efT 

9,q,i 



(T c ) = 0, therefore we add a factor h(T) in the effective densities: 



nf (T) = h(T)n a 



(24) 



with 



h(T) = 1 — exp 



T-T c 
A 



(25) 



To obtain the parameter A = 60 MeV in h(T), we have calculated the pressure correspond- 
ing to these effective densities and compared it with the pressure from lattice simulation in 



Ref. 



27|]. The grand canonical partition function from quarks, antiquarks and gluons has 



the form: 



InZ 



grand 



VTrnj (m a 



(26) 



Here, D a are the corresponding degeneracy factors, D q = Dq = N c Nf(2s q + 1) and D g = 
(N% — l)N p which contain N p = 3 as the number of polarizations of a massive gluon, and 
s q = 1/2 as the spin of the quark. 
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Figure 9: The calculated net entropy S(T) at Nf = 2 (full drawn curve) of a QQ pair at large 
separation is shown as a function of T/T c with T c = 200 MeV. The entropy of the plasma without 
the Q and Q pair is subtracted. Our calculation is compared with a function interpolating the 
lattice data [3, (dashed curve) for better visibility. 



The pressure of free particles 



T 



N 



2miK 2 ( ^ ) + ^miK 2 ( ^ ) + ^±m\K 2 I ^ ) ) (27) 



m 



N 



m. 



6T* 



would not vanish below T c . To have a vanishing pressure at T = T c , we also modify it by 
the confinement factor h(T). Massless quarks and gluons have the density n oc T 3 and the 
pressure P oc T 4 . Therefore, P e ff(T) is approximately suppressed by the factor h(T)^ 3 : 

PMT) = KT f> . ^ ( 2mlK2 (f\ + ^1 <K2 @t\ + ^ m ? A - 2 (^)) . (28) 



In Fig. [3 we plot the ratio P{T)/T A from Ref. [27| and the ratio P cfr (T)/T 4 from Eq. (|28j). 

In the second-order free energy P 2 , Eq. ( f22l) . we include the effective densities n„ ff , 
Eq. (1241) , modified by /i(T). In Fig. [9] and Fig. dOj we have plotted the entropy S and the 
internal energy U following from the standard thermodynamic relations S(T) = —3F 2 /dT 
and U(T) = F 2 (T) +TS(T) for two flavors. For the critical temperature we choose T c =200 
MeV for two flavors [ig|. 

As one can see, S(T) and U(T) above T c are much smaller than in Refs. [si, [jj]]. We 
observe a difference AU — —0.5 GeV for temperatures larger than 1.1T C between our values 



17 



1 1 1 1 r 

our 2 flav. calc. 

interpolation of lattice data 

2 flav. lattice QCD 




Figure 10: The calculated net internal energy U(T) in GeV (full drawn curve) of a QQ pair at large 
separation is shown as a function of T/T c with T c = 200 MeV. The internal energy of the plasma 
without the Q and Q pair is subtracted. Our calculation is compared to a function interpolating 
the lattice data [<], (dashed curve) for better visibility. 

of the internal energy and those of the lattice data, similarly to the result of the unquenched 
calculation at T < T c . The trend from the lattice simulations that S(T) and U(T) decrease 
with T, is also seen here. For the case of three quark flavors, we show the results in Figs.fTTl 
fT2l The calculation of the entropy for three flavors is close to the one for two flavors. A sim- 
ilar result is obtained in the lattice simulations. However, there is one qualitative difference: 
Near T c two flavor lattice simulations give a larger S(T) and U(T) than simulations with 
three flavors. Our calculation yields opposite results. They arise from the flavor dependence 
of the densities and the Debye mass, which partially cancel and lead to an approximate 
independence of S(T) and U(T) on the number of flavors. 

We can enhance the internal energy and entropy by modifying our system in such a way 
that all quarks and antiquarks are bound in color-octet states and interact as such with the 

n 

QQ-p&ir [28fl. This yields an estimate for the influence of bound states on our calculations. 

For simplicity we assume that the qq bound states have a mass of 2m q . The corresponding 
degeneracy factor is now (N% — l)NjN s = 16Nj (N s = 2 for spin or 1). Therefore, the 
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Figure 11: The calculated net entropy S(T) at Nf = 3 (full drawn curve) of a QQ pair at large 
separation is shown as a function of T/T c with T c = 193 MeV. The entropy of the plasma without 
the Q and Q pair is subtracted. Our calculation is compared with a function interpolating the 
lattice data [3, (dashed curve) for better visibility. 



density is given as (cf. Eq. (T23l) ) 



n, 



^9 00 1 /9m 

—TNjm q -K 2 



71 



n=l 



11 



T 



11 



(29) 



As before, we obtain the effective density through the multiplication by h(T), where the 
constant A is again fitted to reproduce the pressure, similarly to Eq. f|28l) . This yields 
A = 80 MeV for two flavors. The color factor c q q = 1/2 is the same as for gluons, and one 
can again determine entropy and internal energy from Eq. ( l2~TI ). These quantities are shown 
in Figs. E3] and H 

In the vicinity of T c the entropy remains smaller than in the lattice calculations, but for 
higher temperatures it is in good agreement with the lattice data. In contrast, the internal 
energy is not very much different from the one without qq bound states. This shows that our 
calculation of the internal energy in second order perturbation theory is not changed very 



much by the presence of bound states in contrast to the otherwise 
effects of bound states in the quark-gluon plasma. 



28, 2t 



important 
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Figure 12: The calculated net internal energy U(T) in GeV at Nf = 3 (full drawn curve) of a QQ 
pair at large separation is shown as a function of T/T c with T c = 193 MeV. The internal energy 
of the plasma without the Q and Q pair is subtracted. Our calculation is compared to a function 
interpolating the lattice data 0, lj| (dashed curve) for better visibility. 



V. SUMMARY 



In the present paper, we have analytically addressed the thermodynamics of a static 
quark-antiquark pair at large separation R > 1.5 fm in the vicinity of the deconfine- 
ment phase transition. In quenched QCD for T < T c , the quark-antiquark string passes 
through valence gluons and forms a gluon chain. We have derived the resulting temperature- 
dependent string tension, Eq. ([6]), which vanishes at the critical point as \JT C — T. In un- 
quenched QCD, for T < T c heavy-light mesons and heavy-light-light baryons are produced 
due to string breaking. The internal energy and entropy of our calculation fit quite well the 
corresponding lattice data. Finally, we have calculated the interaction energy for T > T c in 
unquenched QCD of the Debye-screened static quark-antiquark pair with the constituents 
of the plasma in second order thermodynamic perturbation theory. The change of the free 
energy depends on the effective density of quarks and gluons in the plasma, which has been 
constrained to vanish at T = T c and is adjusted in such a way that the pressure fits the one 
in lattice QCD for Nf = 2. The entropy and internal energy of the model calculation differ 
from the lattice data. The large entropy near T c cannot be reproduced. The internal ener- 
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Figure 13: The calculated net entropy S(T) (full drawn curve) of a QQ pair at large separation 
interacting with qq bound states and gluons is shown as a function of T/T c with T c = 200 MeV. 
The entropy of the plasma without the Q and Q pair is subtracted. Our calculation is compared 
to a function interpolating the lattice data 0, lj| (dashed curve) for better visibility. 



gies U{T) calculated for the heavy QQ pair are in both cases T <T C and T > T c by about 
0.5-0.6 GeV lower than in the lattice simulations. The trend of a falling S(T) and U(T) with 
increasing T is also seen in our model calculation. The heavy QQ-paiv, immersed at large 
separation in the quark-gluon plasma, becomes less and less relevant at higher temperatures. 
Finally, the influence of possible qq bound states on our calculations has been estimated. 
Qualitatively, they do not change our results very much. If one wants to answer the question 
whether the hot QCD system is strongly interacting, our answer is twofold: Below T c the 
gluonic string and the excited hadronic states definitely represent strongly interacting com- 
posite systems. Above T c we have the screened color Coulomb-potential with its strength 
enhanced by the running QCD-coupling. The resulting calulated entropy, however, cannot 
fully reproduce the entropy from lattice simulations. At large temperatures the agreement 
becomes better. Surprisingly the internal energy of the simplified calculation and the lattice 
simulations do not converge. This divergence should be investigated more thoroughly. 
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Figure 14: The calculated net internal energy U(T) in GeV (full drawn curve) of a QQ pair at 
large separation interacting with qq bound states and gluons is shown as a function of T/T c with 
T c = 200 MeV. The internal energy of the plasma without the Q and Q pair is subtracted. Our 



calculation is compared to a function interpolating the lattice data 
visibility. 



(dashed curve) for better 
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Appendix A: COMPUTATION OF THE SUM FOR THE THREE-DIMENSIONAL 

HARMONIC OSCILLATOR 



The sum over n is computed by means of the formula 



1 d 



1 



oo 1 d 00 

^nexp(-A/3um) = -— — £ exp(-A/3um) = , , 

^ (3uj dX ^ pudAl- exp(-A/3a;J 



(Al) 



for the geometric series. Therefore, every factor of n in the sum can be replaced by — 



and we have 



£(^ + l)(n + l)e-^ 

n=0 V Z 



Id 2 3d \ 00 
i 27^dP~2^dA + 1 J5 e 



\f3ujn 
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Figure 15: Interaction of q and Q in second order 



(1 - exp(-A/9w)) 



3 ' 



(A2) 



Appendix B: COLOR FACTORS 



We shall discuss the color structure of the interaction term and the corresponding color 
factors here. As mentioned in the main text, the symbol C a ® carries the color indices of the 
interacting plasma constituent a and both Q and Q. In the first order, the diagrams are 
proportional to TriC a< ^ (or TriC aQ , respectively), where C a ® is a product of the generators 
in the representations of the constituent a and Q. For example, C qQ = t^t^Sik. Now, the 
effect of Tri is that first order diagrams are proportional to the trace of the generator in the 
corresponding representation of the interacting constituent a, which vanishes. Therefore, 
first-order diagrams do not contribute. 

For the relevant second order, we have squares of these color structures. Let us discuss 
this in more detail for one example and find the corresponding color factor. Therefore we 
consider a quark q from the plasma, interacting two times with Q and assume one-gluon- 
exchange (see Fig. [T5l) . 

The color structure is 



^ip Ir mo^pj rk on 



L om L pi °lr i no i jp u rki 



(Bl) 



where the indices k, I and p correspond to the Q which does not participate in this diagram. 
Furthermore, t a are the usual SU(3) generators in fundamental representation. With the 
use of t^t%i = l/2(5a5jk — 1/35^(5^) we can simplify this to 

1 



-5, 



5 1 

~OijO mn ~^jm^in 



(B2) 
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Figure 16: Interaction of a gluon with Q in second order 
As we want to discuss only the QQ singlet state here, we have to apply the projection 



operator 



3ii 



{P?%i,ik = l^ik - ItU-tfkf- (B3) 

This yields the result 2/9 5 mn . A factor of 3, coming from the trace in Eq. ( fTSl) . is taken into 
account for the densities ( l23l . Therefore it should not be counted twice. There are other 
diagrams in second order, namely interaction of q with Q, q with Q and q with Q which give 
the same result of 2/9. However, there are also diagrams, where a quark q interacts with 
both Q and Q or an antiquark interacts with both Q and Q. These diagrams have an odd 
number of interactions in the 3 representation and therefore receive a relative minus sign. 

For the gluons basically the same happens. There is one diagram containing two inter- 
actions of a gluon with Q and the same for Q. In addition, there are two diagrams for 
interactions both with Q and Q which have a relative minus sign for the same reasons as 
above. As gluons are in adjoint representation, we have to use the corresponding generators. 
Let us again consider an example, where a gluon interacts twice with Q (see Fig. [16]). The 
color structure reads 

ri9Q rigQ _ + b Qx -fbac.dQ^ ■ fdce (VIA') 
^irkpac^rjplce ~ L ri u kp>'J L j r u pl l 'J ■ l ±J ^7 

With the immediate use of the projection operator for the QQ singlet state and the Casimir 
operator in adjoint representation f abc f ebc = 3<5 ae , one finally gets l/25 ae . Again, the trace 
over color indices of the gluon was taken into account for the densities and should not be 
counted twice. 

Diagrams which contain interactions of two plasma constituents do not contribute in 
second order because they are products of first-order diagrams which vanish. 
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Appendix C: CALCULATION OF THE SECOND ORDER F 2 



We now calculate the relevant second order for the free energy 



F 2 = - — (V 2 (v 1 ,...,v N )). 



(CI) 



As we discussed before, only terms in the square which come from one plasma constituent 
a contribute. This yields 



1 



N 



(V 2 ( ri ,..., rjv )) = -]Tq dh, 
v i=i J 



2a 



2 ( exp(-m D |r^ - x Q |) ' 
|ri-xo| 



-2a. 



exp(-m D |r^ - x Q |) exp(-m jD |r i - Xq|) 



-a 



(C2) 



Here, V is a three-dimensional volume (coming from Fq) which will be absorbed in the 
densities. While the first term comes from the Yukawa potential squared, the second term 
corresponds to the mixed term of the square. With the use of Fourier transformation, we 
can proceed to (dj = r — Xj) 



d 3 r, 



(2tt) 3 (2 



{ql + m 2 D ) 2 



c? 3 qi d 3 q 2 ( . exp^qudg) exp(iq 2 dg) 

3(^)S a -(^ r ) ( ^ + m 2 ))2 «.(ft. r ) 

= -- I d 3 qa s (q,T) 2 



-2(4tt) 2 / 



,exp(iq(x Q -Xq)) 



(q 2 + m 2 D ) 2 



(C3) 



for the second term and 



2 1 (PrJa.dri-xgl.T) 

7T J 



exp(-m jD |r i - xq| 

l r * - x qI 



(C4) 



"(g 2 + m^) 2 

for the first one. The sum becomes trivial for each species of plasma constituents and can 
be written as E a c a Nf- With Nf /V = rf^ one can then introduce the densities here. 
Altogether we have 



Oi s (q,Tf 



(q 2 + m 2 D ) 2 



(» 



^q(XQ-Xg) 



)■ 



(C5) 
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Appendix D: COMPUTATION OF THE DENSITY OF FREE RELATIVISTIC 

PARTICLES 

We want to calculate the density of free relativistic particles for quarks (fermions) and 
gluons (bosons). It has the form 

/d?p 1 
TTTs (qj? \ Z i ( D1 ) 
(2-k) 6 exp(pE ) ± 1 



where a stands for q, q or g, E = yp 2 + m 2 , and D a is the degeneracy factor. Explicitly, 
this factor is (iV 2 — l)N p — 8 • 3 = 24 (due to color and polarization) for (massive) gluons 
and N c Nf(2s q + 1) — 3 • Nf • 2 = 6iV/ (due to color, flavors and spin) for quarks. 

Let us begin with quarks. The denominator can be rewritten by means of the geometric 
series: 

I I I e _^ g ( _ 1)ne _^ £o _ (D2) 

n=0 



exp(/3E ) + 1 exp(/?£ ) 1 + exp(-(3E 
Next, we can use the general formula 



mR 

c 

4a — 

o \/7ia m 



00 da ^_ m 2 a _!? e 



(D3) 



to obtain 

l2N f f)(-i)«+i |^e-^ 



n=l 



The integral over p is now Gaussian and gives f^p) • With the use of the formula 

£° dti*- 1 exp ^-t* - ^ = 2 ^ K a {2^6), Re 7 > 0, Re 5 > 0, (D5) 



we get the final result 



«« = 3^gfc& 2 (^). (D6) 



The same result will appear for antiquarks, with m 5 instead of m q . For gluons (bosons) 
with the thermal mass m g the computation is nearly the same. Due to the minus sign in 
the denominator of ( 1D1I ). the factor (—1)" does not appear (and, of course, the factor of 
degeneracy is changed). We have 
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